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Abstract 

We consider a simple renormalizable dark matter model consisting of two real 
scalars with a mass splitting 5, interacting with the SM particles through the Higgs 
portal. We find a viable parameter space respecting all the bounds imposed by 
invisible Higgs decay experiments at the LHC, the direct detection experiments by 
XENONIOO and LUX and the dark matter relic abundance provided by WMAP 
and Planck. Despite the singlet scalar dark matter model that is fragile against the 
future direct detection experiments, the scalar split model introduced here survives 
such forthcoming bounds. We emphasize on the role of the co-annihilation processes 
and the mixing effects in this feature. For ?«dm ~ 63 GeV in this model we can 
explain as well the observed gamma-ray excess in the analyses of the Fermi-LAT 
data at Galactic latitudes 2° < |6| < 20° and Galactic longitudes |£| < 20°. 
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1 Introduction 


Although there is no doubt on the existence of dark matter (DM) which is forming 
about 26 percent of the matter content of the Universe m (see e.g. reviews em 
its fundamental interaction with ordinary matter of the Standard Model (SM) of 
particle physics is a tremendous mystery in physics today. There is however, a 
natural explanation for the present value of DM relic density in terms of the thermal 
freeze-out mechanism of weakly interacting massive particles (WIMPs). Exploiting 
the WIMP paradigm, a large number of theories beyond the SM is developed with 
a DM candidate as a WIMP, we name for instance supersymmetric models with 
R-parity , models with universal extra dimensions, as well as models with minimal 
extension of the SM which is of our interest in this article ehhi. 

All these models can receive stringent constraints on the DM annihilation cross 
section from Planck Jl| and WMAP (2], precise measurements of the DM relic den¬ 
sity, and on the WIMP-nucleon scattering cross section from dark matter experiments 
such as LUX [T3] and XENONIOO [T4l . Moreover, in the case of DM production in 
particle collider experiments, there are measurements such as invisible Higgs decay 
and missing energy-momentum that can put further restrictions on the model pa¬ 
rameter space (151IT8] . 

The new bounds by the coming direct detection experiments such as XENON1T 
which is going to start data collection already in this autumn, will certainly exclude 
many of the current WIMP models. The popular singlet scalar dark matter model 
as the most minimal extension of the SM has been investigated elaborately from 
different points of view in the literature (see for instance [5ll6| 119ll20] ). Although this 
model is fairly successful in various aspects, it is quite in danger to be excluded for a 
wide range of DM mass due to the direct detection experiments that will put stringent 
bounds in the near future [2TJ. If the direct detection experiments are taken seriously 
and one is still interested in the scalar extension of the SM, the next minimal model 
that comes to mind is the two real scalar extension dubbed here under the name of 
scalar split WIMPs. We show in this paper that scalar split model is as good as the 
singlet scalar model with drastically improved features in the direct detection part. 

On the other hand, in the light of the recent confirmed observation of the Fermi- 
LAT extended gamma ray excess, many investigations have directed towards possible 
explanation of the gamma excess. Assuming that the galactic gamma excess produced 
as a result of DM annihilation in the galactic center, it is then found in a number 
of models that DM annihilation cross section of order rs_/ 10 26 cm 3 s 1 with DM 
mass in the range 30 — 50 GeV can explain the excess, see as examples |22H49| . and 
see [50H53] for scenarios with lighter DM. Later it was found in [55H5lj that DM mass 
of ~ 35 — 165 GeV decaying into b quark pair and DM mass large enough to decay 
into W + W~, ZZ, hh, tt pairs can be fitted satisfactorily to the Fermi-LAT data. 

In this paper we consider a minimal extension of the SM with two additional real 
scalars denoted by S\ and S 2 , which are SM gauge singlets and interact with the SM 
particles via a Higgs portal respecting the Z 2 symmetry under which the new scalars 
are odd and all the SM particles are even. This model suggests two scalar WIMPs 
with a mass splitting 6 where only the lighter component is stable and the heaver one 
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is an unstable state jS]. The viable parameter space constrained by the limits from 
the observed DM relic abundance, direct detection bounds as well as invisible Higgs 
decay width is studied in this work. We also show that it is possible to find regions in 
the viable parameter space which can explain the galactic gamma ray excess observed 
by Fermi-LAT. 

The rest of the paper has the following structure. In section [2] the scalar split 
model is introduced and the relevant free parameters are discussed. In section [3] we 
discuss how to calculate the relic density using the DM annihilation cross sections. 
Section [4] is devoted to calculations on the Higgs decay to two WIMPs and invisible 
Higgs decay width is provided in terms of the mass range of the DM candidate. 
Moreover, the viable parameter space constrained by the DM relic density observation 
as well as invisible Higgs decay width are studied. Elastic scattering cross section 
of DM-nucleon is computed as a function of DM mass in section [5j taking into 
account the limits from relic density observation and direct detection experiments. 
In section [6] we find how it is possible to explain the inner galactic gamma ray excess 
within the constrained model parameters. We finish in section [7] with the conclusion. 


2 Scalar Split WIMPs 

We consider a renormalizable extension to the SM with two new real scalar fields 
denoted by S± and S 2 . These new fields may have small mass splitting and transform 
under 7L 2 symmetry as Si —» —Si. The full Lagrangian consists of 

£ = £sm + £dm + -^int • (i) 


The Lagrangian £dm incorporates only the WIMPs particles as 


1 


1 


£dm = t;{9u,S 1 ) + -(9 m 52) — 
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'S 2 - 


m. 
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2 

2 C 2 
d 2 ~ 


A'i q4 A4 q4 
T ^- T ^2 


( 2 ) 


We could in principle have the interaction term A s S 2 S 2 in the lagrangian ([2]). We 
will turn to this point later in this section. 

In addition, respecting the Z 2 symmetry, WIMPs interaction with the SM parti¬ 
cles are considered through a Higgs portal such that 


£ int (Si ,S 2 , H) = Ai SfH'H + \ 2 S 2 H^H + A 12 Si S 2 H*H . (3) 


The SM-Higgs potential is also given by 

V H = » 2 h H'H + \ h (H'H) 2 . (4) 


The Higgs field is a SM SU(2)l scalar doublet which develops a non-zero vacuum ex¬ 
pectation value (vev) which results in the electroweak spontaneous symmetry break¬ 
ing. We then parameterize H as 



v + h ) ’ 


( 5 ) 
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where v = 246 GeV. 

We can choose a basis in which (Si) = (S 2 ) = 0. The minimization conditions of 
the total potential are 


dV 

dV 

dV 

dH 

(H)=v/V2 dSi 

<Si)=o dS 2 


(S 2 )=o 


= 0. 


( 6 ) 


These conditions provide us with some relations between the parameters. We work 
them out and identify the entries of the mass matrix. From condition ^ | ^ = 0 we 
get the relation 


9?h ~ ■ 


(7) 


From the other two minimization conditions we get no more relation. We also get 
the following results for the entries of the mass matrix 


and 


2 d 2 V 2 . 2 2 d 2 V 2 x 2 

m S l = 332 = m l + A W > m S 2 = -gg2 = m 2 + ^2V , 


m Si,S 2 


d 2 V 

dS 1 dS 2 



( 8 ) 

(9) 


We then indicate the two fields H 1 and H 2 as the mass eigenstates by introducing 
the mass mixing angle 9 


H 1 = sin 9 Si + cos 9 S 2 , 

Hi = cos# Si — sin# S 2 , (10) 


where, 


tan 9 = 


with y = 


2m 


Si,S 2 


2 2 
171 S 2 ~ m S\ 


1 + y/l +y 2 

The two neutral scalars H\ and H 2 have the corresponding mass eigenvalues as 


( 11 ) 


m H u H 2 


m li + "4 ± "4 - m li 


( 12 ) 


We assume that > mH 2 an d therefore H 2 is the stable DM candidate. It is then 
possible to obtain relations for m\ and m 2 in terms of physical masses and couplings 

777-1 = m 2 Hl sin 2 9 + m 2 Hi cos 2 9 — X\v 2 , 

m 2 = m 2 Hl cos 2 8 + m 2 H2 sin 2 9 — \ 2 v 2 . (13) 

Moreover, from eq. (j9j) , eq. m and eq. (fT2|) we can express the coupling A 12 in terms 
of the masses m h 1 and mn 2 and the mixing angle 9 , 


A 12 



m 2 H2 ) sin 29 . 


(14) 
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Figure 1: Three body decay of the scalar Hi into the scalar DM and a fermion pair. 

We now turn to the point we made after equation (j2j). If we rewrite the lagrangian 
Q in the basis of mass eigenstates H\ and H 2 , then it includes a term for interacting 
DM with its parter as 6 sin 2 6 cos 2 Taking rrijj l > m # 2 it means that 

the co-annihilation process H 2 H 2 —» H\H\ kinematically is not allowed. Adding the 
AsS'xS ' 2 term do not introduce any new interactions in the lagrangian after going to 
the mass eigenstate basis, although it modifies the strengths by the new coupling A s . 
For instance, the term above changes as (6 sin 2 9 cos 2 0{ A 3 + A 4 — 12A S ) + 6 A S ) HfH%, 
which again do not contribute in the relic density computation. Therefore, the term 
S' 2 S'! merely enlarges the dimension of the parameter space by one. To stay in the 
most minimal scenario possible we assume that A s = 0 in this paper. 

We therefore can take in the present model seven independent parameters as m/q, 
m# 2 , Ai, A 2 , A 3 , A 4 and 9, while the coupling A 12 is then fixed by the relations in 
eq. CD and eq. m- The vacuum stability of the total potential restricts the model 
parameters. In this regards, we find at tree level the bounds 


771 2 + 77T , 2 "p (Al + A 2 )u 2 > 0 , 
+ (mfXi + rn^^v 2 + A1A2U 4 > Xf 2 v 4 ■ 


(15) 


In addition, the perturbativity of the model requires the upper bounds on the cou¬ 
plings, | Ai| < 47 t. 

When the small mass splitting is the case then the heavy component WIMP can 
decay into an off-shell Higgs and the light partner as H\ —>• H 2 I 1 where h itself decays 
successively into a fermion pair as h —>• //. The Feynman diagram for the decay is 
shown in Fig. [D 

It is necessary to have an estimate on the life time of the heavy component over 
the restricted parameter space to know whether or not it has any contribution on the 
DM relic abundance. We provide here the formula of the double differential partial 
decay width for Hi{k) H 2 {p 2 .) fipi) /(P 2 ) 


d?T 


3m 2 [(Ai — A 2 ) sin 29 + A 12 cos 26 ] 2 r t 



- m H 2 ~ 4m /l 


dt du 


1287 T 3 rrifl i 


where the mandelstam variables are t = (pi + P2) 2 and u = (p2 + Pz) 2 ■ 
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Figure 2: The Feynman diagrams for the DM (co)-annihilation into SM hnal states. Dia¬ 
grams with more than two particles in the final state are not shown. 

3 Dark Matter Relic Abundance 

Assuming that DM particles have been in thermal equilibrium in the early Universe, 
the present density of DM depends somehow on the so-called freeze-out temperature, 

Tf , the epoch in which dark particles become non-relativistic and go out of the 
equilibrium. At freeze-out temperature the annihilation rate of DM falls off below 
the Hubble expansion rate. On the other side, due to the low budget of the kinetic 
energy, the DM production reactions get suppressed. The relic density of DM is 
computed by solving the Boltzmann equation(s) for the time evolution of DM number 
density, udm- In the model under consideration, there are two new scalars beside the 
SM particles that their number density evolutions are relevant in order to obtain the 
DM relic abundance. We assume that H 2 is the lighter component and thus is stable. 

We therefore consider H 2 as our DM candidate with mass m# 2 that rriH 1 > ■ So 

the heavier scalar H\ can undergo the decay H\ —>• H 2 + SM. 

Annihilation reactions are one type of processes that change the number density 
(ni and 712 ) of our species here. The possible annihilations of H\ and H 2 to SM 
particles are depicted in Fig. [2j As it is evident from the Feynman diagrams, anni¬ 
hilation reactions into SM fermion pairs, W + W~ and ZZ occur via s-channel while 
annihilation into SM-Higgs pair is possible through s-, t- and u-channel. An anni¬ 
hilation process in which DM particle annihilates together with H\ is the so-called 
co-annihilation reaction. Another type of reaction that changes the number density 
is the decay process of the heavier component, i.e., H\. In principle, the abundance 
of H 1 and FD are determined by solving two coupled Boltzmann equations. The two 
Boltzmann equations can be written in a single Boltzmann equation with an effective 
(co-)annihilation cross section [58h60| . 

^ = -3 Hn - (a e ffv) ( n 2 - n 2 eq ) , (17) 
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Figure 3: The dependency of the relic density on the mass splitting 5 for a wide range of 
DM mass. 

where n = nn 1 + n # 2 and 

a eS = — ( cr 2 2 + On fl + e~ 2S/T + 2 cti 2 (l + , (18) 

deft \ V m H 2 J V m H 2 J J 

where <t 2 2 , cth and 012 stand for (co-)annihilation processes H 2 H 2 —> SM SM, H\H\ —> 
SMSM and H 2 H\ —> SMSM respectively with g e g- = 1 + fl H—— ) ^ e~ 5 ^ T . 

The expression (n e fjn) indicates thermal average over effective annihilation cross 
section x relative velocity at temperature T. In appendix A we present the formulas 
for annihilation cross sections of dark matter candidate in four possible channels. To 
confirm our analytical formula we employ the program CalcHEP }61] which in turn 
requires implementation of our model into the program LanHEP [62). To perform 
the analysis for the DM relic abundance we need to solve numerically the Boltzmann 
equation. To this end, we utilize the program MicrOMEGAs |63j for our model. 

As explained earlier we have two choices for a set of independent parameters 
we would like to place the constraints on. Notice that the couplings A 3 and A 4 
do not show up in DM annihilation cross sections at the tree level, however these 
couplings appear through the strength of the vertex HfH^- We fix the two couplings 
as A 3 = A 4 = 0. Therefore, one possibility is choosing the set of parameters {rriHi, 
mn 2 , Ai, A 2 , A 12 } and the other option is the set {mn 1 , i^h 2 , Ai, A 2 , 0}. In our 
analysis we choose the second set and apply the relation in eq. m to obtain the 
coupling A 12 by fixing the mixing angle 6. 

Let us define the mass splitting as <5 = A?ni 2 = m f/, — mn 2 ■ Taking into account 
the co-annihilation processes, we check numerically the dependency of the relic den- 
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sity on the mass splitting 5. For a point in the parameter space with Ai = 0.56, 
A 2 = 0.33 and sin# = 0.1, the results are compared in Fig. [3] for <5 = 1 GeV, 4 GeV, 
40 GeV and 100 GeV. Since the co-annihilation effects are larger for smaller value of 
5, for the present model, it is evident from Fig J3]that the relic density is reduced by 
the co-annihilation effects. 


4 Invisible Higgs Decay 


The DM candidate in the scalar split model interacts with the SM particles via SM- 
Higgs mediator. It also opens up the possibility for the 125 GeV Higgs to decay 
into the new scalars. Constraints on the model parameters are placed by requiring 
the invisible Higgs decay to be consistent with the Large Hadron Collider (LHC) 
measurements. The total decay width of 125 GeV Higgs decaying into SM particles 
is ~ 4.1 MeV [64] which get enhanced by three invisible decay width of the SM-Higgs, 
h —» H\. h —» H 1 H 2 and h —>• H 2 H 2 . Given an experimental upper limit for the 
invisible branching ratio for the Higgs boson as Ti llv /(Ti nv + T v ) ~ 0.35 [18] we put 
a bound on the total invisible decay width as ri°* al < 2.15 MeV. On the other hand, 
the total invisible decay width in this model is saturated by three possible decays of 
the Higgs: 


r &(h Hi 


r f 2 v (h h 2 h 2 ) = 


8tt nrih 


1 (Ai cos 2 9 + X 2 sin 2 9 — A 12 sin 9 cos 9) 2 4 m 2 H2 


8itmh 


4m ^i yn 

(19) 

m 2 h 


4 rr 4 2) 1/2 

( 20 ) 


m 


and 


TZ(h = 


v 2 [(Ai — A 2 ) sin 29 + A 12 cos 2 0y 
8 nml 


[ml ~ (■ m 2 Hl + m 2 H2 ) 2 } 1/2 [m 2 h - (m 2 Hl - m^ 2 ) 2 ] 1/2 . (21) 


The invisible Higgs decay width depends on DM mass m # 2 and 5 as the following: 


-ptotal 
^ inv 

p22 

^ inv 

when 

m h 

2 

5 nih 

~ 2 <mH * < ~2~ 

ptotal 
^ inv 

_ p22 1 pl2 

inv ' 1 inv 

when 

m h 

2 

-5 < m H2 < — 

ptotal 
^ inv 

_ p22 1 pl2 1 pll 

inv ' inv ' inv 

when 


^ m h 

m H 2 < ~y ~ ' 


( 22 ) 


Let us now begin with our probe over the parameter space of the model. To 
proceed we put together the constraints imposed on the parameter space from relic 
density analysis and the invisible Higgs decay width. This has been done for two 
values of the mass splitting 5 = 1 GeV and S = 100 GeV in order to investigate the 
role of this parameter on the viable space and confronting that with the singlet scalar 
dark matter case. 
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Figure 4: Shown are the allowed DM mass in the viable parameter space respecting the 
relic density and invisible Higgs decay width constraints for left ) 5 = 1 GeV right) S = 100 
GeV. 


We have generated random values for DM mass with 40 GeV < m// 2 < rrih/ 2, 
—1 < Ai < 1, and taking A 2 = Ai/5 and sin 0 = 0.1. Using the combined results from 
WMAP and Planck for the present DM relic density, the results exhibited in Fig. [J] 
show the viable parameter space for different values of the mass splitting parameter 
5. The region colored in blue is excluded by the invisible Higgs decay width. There 
are some comments in order for the Fig. [U 

First, for DM mass below m^/2, only DM annihilation into fermions mediated 
via SM-Higgs are potentially allowed processes, thus one expects enhancement on 
the cross section near the SM-Higgs mass resonance. It is evident from Fig. [4] that 
the resonance occurs around DM mass mjj 2 ~ m^/2 ~ 62 GeV where the coupling 
A 2 (as well as the annihilation cross section) takes its smallest value, while A 2 grows 
up for DM masses smaller or greater than the resonance mass 62 GeV. Moreover, it 
is seen that for both mass splittings 5 = 1 GeV and 5 = 100 GeV all DM masses 
smaller than mj 2 ~ 56 GeV are excluded. Finally, the range of the allowed DM mass 
in the scalar split model is almost the same as that of the single scalar dark matter 
model where 7 tidm < 55 GeV is excluded by the LHC bounds on the invisible Higgs 
decay width [21] . 


5 Direct Detection 

Direct detection experiments are designed to study the unknown nature of DM inter¬ 
action with ordinary matter. In these experiments the attempt is to measure the en¬ 
ticing event rate for the DM scattering off the target nuclei in the detector. Although 
the present results from DM experiments such as LUX [[13] and XENON100 m show 
no evidence for DM interactions, they offer an impressive upper bound on the spin- 
independent DM-nucleon elastic scattering cross section. We will apply these findings 
in the following to constrain further the parameter space of our model which is al¬ 
ready restricted by the limits from WMAP and Planck. To this end, we need to 
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Figure 5: The relevant Feynman diagram for the WIMP-nucleon elastic scattering. 

calculate the elastic scattering of WIMP-nucleon. In the present particular model 
the interaction of DM with nucleon occurs through a fundamental interaction of DM 
with quark which is mediated by the SM-Higgs, where the relevant Feynman dia¬ 
gram is depicted in Fig. (5[ The effective Lagrangian responsible for the DM-quark 
interaction is, 


£ eff = a q H 2 H 2 qq , 


(23) 


where, the coupling a q is given by 
im 

a q = —§■(Ai cos 2 6 + \ 2 sin 2 6 — \\ 2 sin 6 cos 9 ). (24) 

m h 

To find the elastic scattering cross section we can invoke the assumption that in the 
limit of vanishing momentum transfer it is possible to replace the nucleonic matrix 
element including quark current with that containing nucleon current up to some 
proportionality factor [651168] . see also [69] • We arrive at the final result for the 
spin-independent (SI) cross section of DM-nucleon as 


N a N^N 
°SI — 


-Km 


DM 


(25) 


where /tat is the reduced mass of the DM-nucleon system and the factor a n depends 
on the scalar couplings f^ q and as 


q;jv = m N 


E 

q=u,d,s 



, 1 f N \ ^ a <l 

+ 27 at 9 2 - , m 

q=c,bit * 


(26) 


In our numerical calculations we use the following values for the scalar couplings 


fP = 0.0153, fP = 0.0191, fP = 0.0447. (27) 

We compute the elastic scattering cross section for about 10 6 points in the pa¬ 
rameter space by random generation of the relevant parameters with — 1 < Ai < 1, 
—1 < A 2 < 1, 40 GeV < j«dm < 3 TeV and sin# = 0.1. We report on our numerical 
results for the elastic scattering cross section of DM-proton in Fig. [6] and Fig. [7] con¬ 
sidering four different mass splittings in the model, namely 5 = 1,10,40,100 GeV. We 
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Figure 6: Spin-independent DM-nucleon elastic scattering cross section are shown as a 
function of the DM mass and comparison has made with the latest results from LUX and 
XENONIOO experiments and the future experiment XENON1T. In the left panel the mass 
splitting is 5 = 1 GeV and in the right panel 5 = 10 GeV. The vertical color spectrum 
indicates the size of A 2 . We have chosen for the couplings as — 1 < Ai < 1 and — 1 < A 2 < 1. 




Figure 7: The same as in Fig. [6] with the left panel for the mass splitting § = 40 GeV and 
the right panel for 5 = 100 GeV. 
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Figure 8: The elastic scattering cross section with sin (9 = 0 for the parameter space 
bounded by the present Universe DM relic density. 

have exploited the viable parameter space fulfilling already the relic density bound 
to obtain numerically the elastic scattering cross sections for a wide range of DM 
mass. These results have been plotted against the experimental bounds provided by 
LUX and XENON100 and the estimated bound for the future XENON1T. 

In all plots it can be seen easily that the cross section falls off suddenly at the 
Higgs mass resonance region as pointed out before in section |5J It is for the resonance 
mass, i.e., around 62 GeV and for moM > 4 TeV that the singlet scalar model m 
can evade the future direct detection bounds. The success for the scalar split model 
is that it not only can evade the LUX and the XENON100 constraints but there are 
many DM candidates for which the values of the elastic scattering cross sections go 
even much below the future direct detection experiments such as XENON1T for a 
quite wide range of DM mass. This feature is true for all mass splittings 6. However, 
the viable parameter space is slightly sensitive to the mass splitting as seen in Figs. [ 6 ] 
and [71 When 5 is small, i.e., when DM scalar and its partner in the model have more 
or less the same mass, the viable DM mass can be in the range ~ 50 — 200 GeV for 
<5 = 1 GeV if we consider the constraints imposed by XENON1T. With increasing 
the mass splitting 5, the viable DM mass is limited to values in the range ~ 125 — 200 
GeV. It seems therefore that for small enough <5 the parameter space is greater than 
that with much bigger 6. Looking at Fig. (3[ note that for the viable space, i.e., for 
moM ~ 125 — 200 GeV, the relic density does not change considerably going from 
<5 = 1 GeV to (5 = 100 GeV. On the other hand, the value of the splitting mass 5 can 
make a big change in the relic density for ?bdm < 60 GeV and todm > 400 GeV. 

We claim that the improvement in the spin-independent DM-nucleon scattering 
cross section for the present scalar split model of DM stems from two distinct effects. 
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One is the contribution of the co-annihilation processes and the other one is what we 
call it mixing effect. These two effects do not exist in the singlet scalar DM model. 
To clarify this statement, we repeat the computation for the DM-nucleon scattering 
cross section for sin 0 = 0 where the contribution from the co-annihilation processes 
in the relic density and the mixing effect are both absent. The results for ogi that 
respect the relic density bound is plotted in Fig. [HI This figure actually accounts for 
the similar results presented in m for the singlet scalar DM model. It is evident 
from Fig. [8] that the feature in Fig. [6] and [T| disappears in the case sin0 = 0. There is 
a simple explanation why the new feature is not possible in the single scalar model. 
The reason hinges in the fact that in this case, both the annihilation cross section and 
the DM-nucleon scattering cross section are proportional to one common parameter, 
Ai, as can be seen from the formulas provided in the appendix [Aj So that it is not 
possible to get simultaneously a quite small value for the DM-nucleon scattering 
cross section and a large enough annihilation cross section suitable for predicting the 
observed relic density. 

The question we would like to address here is that why the new feature for large 
5 starts appearing in the region with ttidm >125 GeV when 0 / 0. We know that 
for 771 dm < 125 GeV, only the processes, H 2 H 2 —>• / f,W + W ~, ZZ contribute to 
the total annihilation cross section. Looking at the relevant formulas given in the 
appendix El we find out that both the annihilation cross section and DM-nucleon 
scattering cross section are proportional to one common parameter, /?, where (3 = 
Ai cos 2 0 +A 2 sin 2 0 —A 12 sin 0 cos 0. Therefore, we can apply the same line of reasoning 
as we did in the singlet model to understand why in the region with ttidm < 125 
GeV we see the same prediction as the one in the singlet model. Of course, for small 
mass splitting, i.e., (5 = 1 GeV the co-annihilation effects are sizable such that for 
771 dm <125 GeV, the spilt scalar model and singlet model show different predictions 
for the DM-nucleon scattering cross section. Now, when ?hdm gets values larger than 
125 GeV, the process H 2 H 2 —>• hh starts dominating the total annihilation cross 
section. Here we expect the mixing effects show up. In this region it is totally possible 
to find quite small values for f3 and hence small DM-nucleon scattering cross section. 
At the same time having large values for a results in large enough annihilation cross 
section, where a = (Ai — A 2 ) sin 20 + A 12 cos 20. To see this point, we need to look at 
eq. (1321) in the appendix El We see that even when f3 is small, the annihilation cross 
section can be large enough since a is not necessarily small and terms involving a 
will dominate the annihilation cross section. We have justified this latter claim in 
our numerical computations. When DM mass is larger than ~ 188 GeV such that 
the process H 2 H 2 —>• hhh becomes kinematically possible, the aforementioned mixing 
effects discussed above are plausible, however, their strength would depend on the 
size of the couplings Ai, A 2 and <5. Heavier DM with mass 7 tt,dm <1 250 GeV will open 
the new channel H 2 H 2 —>• hhhh and we can see its small effects in Fig [6] and Fig [3 

We redo our computations with couplings in the range —5 < Ai < 5 and —5 < 
A 2 < 5 for <5 = 1,10,40,100 GeV. Our results given in Fig. [9] for <5 = 1,10 GeV 
indicate that for larger values of the couplings Ai and A 2 , DM candidates which can 
evade XENON1T constraints are extended to masses up to ~ 1000 GeV. The DM 
candidates are extended to masses up to ~ 500 GeV in case <5 = 40,100 GeV, as 
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Figure 9: The same as in Fig. [6] with —5 < Ai < 5 and — 5 < A 2 < 5. 



Figure 10: The same as in Fig. dwith — 5 < Ai < 5 and —5 < A 2 < 5. 

depicted in Fig. [TOj Comparing our results in Figs. 191101 with those in Figs. 16171 we 
realize that the effects associated with the processes H 2 H 2 —>• hhh and H 2 H 2 —>• hhhh 
become sizable for |Ai| > 1 and |A 2 1 > 1. 

For the current model, inelastic WIMP-nucleon interaction begins contributing 
for 5 ~ KeV which is far smaller from the limits we have considered in this work i.e., 

5 ~ GeV. 

6 Gamma-ray Emission from DM Annihilation 

The analysis of Fermi Large Area Telescope (Fermi-LAT) data [70] (see m for the 
recent report) triggered by the authors in [72773] and continued by several groups 
lEMd-EE] revealed an excess in the gamma-ray from the center of the Milky 
Way or Galaxy Center (GC), hence dubbed Galactic Center Excess (GCE). The 
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gamma-ray emission produced by the millisecond pulsars in the galaxy center can 
only contribute in 5-10% of the excess observed [82] @. Sources such as cosmic ray 
interactions are disfavored as well [ 83] . 

On the other hand, surprisingly the morphology and the spectrum of the GCE is 
well fitted when the dark matter annihilation into standard model particles is added 
in the background model used in the analyses. All diffuse background models where 
include the WIMP as a component agree in morphology. However, the position 
of the gamma-ray peak in the energy spectrum and the mass of the dark matter 
annihilating into SM particles varies by considering the systematic uncertainties in 
the background model [ 551 56]. 

The dark matter candidate depending on its mass can annihilate into leptons, 
quarks, the gauge bosons and the Higgs boson. The gamma-ray is then produced 
through the cascade decays of these particles to neutral pion by the hadronization 
of the quarks, also through the bremsstrahlung of the charged gauge bosons and 
leptons. Among these processes the gamma-ray from the pion decay is dominated 
compared with the gamma emission from bremsstrahlung. The differential gamma- 
ray flux produced by a single W, Z, the Higgs boson and the top quark is depicted in 
Fig. 1 of [55] • It can be easily seen that the peak of the spectrum is moving towards 
the higher energies for heavier particles. 

It was believed formerly (see e.g. [78lf80] ) that dark matter candidates with masses 
being only in the range of 30 GeV < moM < 50 GeV decaying into bb give an 
acceptable fit with the excess observed in the Fermi data. In the recent works however 
it is argued that taking into account the systematic uncertainties in the analysis of 
the Fermi data not only the mass range of dark matter for bb channel is enlarged into 
35 GeV < tu-dm < 165 GeV but also larger dark matter masses in annihilation to 
WW, ZZ, hh, and tt can be fitted well enough with the data m- Additionally, it 
is pointed out in m that DM mass up to about 74 GeV decaying into b quark pair 
and also DM annihilation into non-relativistic hh can fit well to the Fermi data. 

We show that the gamma-ray excess in our scalar split model can be explained 
well. To this end, we obtain the photon flux produced by dark matter annihilation 
where the allowed values for the couplings are taken from the viable parameter space. 

The gamma-ray flux is determined in terms of the annihilation cross section 
(<7u)ann) the mass of the annihilating dark matter todm, the ganrnra-ray spectrum 
generated per annihilation dN^/dE 7 and the density of dark matter p in the region 
of interest (ROI) is 



(28) 


The density of dark matter in the Milky Way galaxy is assumed to be spherically 
symmetric. The density distribution is then a function of r and is described by the 
generalized Navarro-Frenk-White (NFW) halo profile [84] : 



(29) 


1 Although a recent paper [85] associates GCE to the point-like sources such as millisecond pulsars. 
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Figure 11: Shown is the gamma-ray flux multiplied by energy squared from annihilating 
dark matter computed in the scalar split WIMPs model for dark matter mass ?? t-dm ~ 63 
GeV and S = 100 GeV. The black error bars accompanied with correlated systematic errors 
is the obtained flux from Fermi-LAT data [ 56] , 


where r s = 20 kpc is the scale radius, p & = 0.3 GeV/cm 3 is the local dark matter 
density at r 0 = 8.5 kpc and r is the distance from the center of the galaxy to the point 
where the dark matter annihilation occurs. The parameter 7 is the slope parameter 
being 7 = 1 for the standard NFW. In our calculations we take 7 within the interval 
7 = 1.2 — 1.3 used in the literature. 

We find out that within the parameter space confined by relic density and direct 
detection in our DM model, there can be found regions producing gamma excess 
that are compatible with the fluxes provided by the Fermi data. We have used 
MicrOMEGAs package for computation of the gamma-ray flux in our particular 
model with dark matter mass uidm ~ 63 GeV, mass splitting 6 = 100 GeV, Ai = 
1.17 x 10 -2 , A 2 = 6.07 x 10 _1 and sin 9 = 0 . 1 . 

In Fig. [lT]we present our results for the gamma-ray flux multiplied by the gamma 
energy squared for two slope parameters 7 = 1.25 and 7 = 1.30. In this figure 
it is shown the prediction of the scalar split WIMPs model for the gamma excess 
from annihilating dark matter of mass ~ 63 GeV with the total annihilation cross 
section (<x) ann ~ 4 x 10~ 26 cm 3 s _1 and DM-nucleon elastic scattering cross section 
cr e i ~ 8 x 10 ~ 10 pb, to be compared with the excess observed from the Fermi-LAT 
data. Comparison made with the data analysis provided by [56 ] at high Galactic 
latitudes 2 ° < | 6 | < 20 ° indicates the validity of the current model in explaining the 
Galactic gamma excess. 
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7 Conclusions 


In this paper we have employed a simple model of dark matter called Scalar Split 
WIMPs with two scalars H\ and H 2 interacting with SM particles through the Higgs 
portal. Depending on the mass splitting 6 = rn f/, — mn 2 and the couplings in 
the model, the decay rate of the heavier scalar H\ changes. For the parameter 
space explored in this work, the H\ decay rate is much smaller than the age of 
the Universe. Therefore we have only one scalar H 2 that contributes in the DM 
relic abundance. The model possesses seven free parameters out of which only five 
parameters run 1 , m// 2 , Ai, A 2 and 0 enter into the annihilation and co-annihilation 
cross section computations. The mass splitting parameter <5 may change the viable 
parameter space in various computations. 

We have examined our model with four observational and experimental bounds 
imposed by invisible Higgs decay, the amount of dark matter abundance, the limits 
put on DM-nucleon cross section by direct detection experiments, and the gamma 
excess found by the analyses on the Fermi-LAT data. 

There is an important characteristic for the scalar split WIMP model that dis¬ 
tinguishes that from the vastly studied singlet scalar models. In case we choose 
—1 < Ai < 1 , — 1 < A 2 < 1 , there can be found viable regions in the parameter space 
beyond the resonant region in the scalar split model with ttidm in the range around 
57 — 200 GeV for <5 = 1 GeV in Fig. [6] and around 125 — 200 GeV for <5 = 40,100 
GeV in Figs. [71 which evade the future experiment bounds on the WIMP-nucleon 
elastic scattering cross section. The viable region is extended to higher DM mass up 
to ~ 1000 GeV if we choose —5 < Ai < 5, — 5 < A 2 < 5. For the DM candidates 
in the viable space that respect the XENON1T and LUX bounds, we have inferred 
that the presence of the co-annihilation (when the mass splitting is small) and the 
mixing effect in the scalar split model play a critical role in the new feature so that 
the model can evade easily the bounds from the future direct detection experiments. 

This feature is absent in the singlet scalar model where as shown in [21] the 
viable region in the parameter space is confined in the resonant region with ttt-dm 
in the range 57 — 62 GeV or for ?udm <1 4 TeV by the forthcoming direct detection 
experiments. 

We observe that changing the mass splitting <5 has almost no effect on the invisible 
Higgs decay width. Furthermore, the behavior of the relic density for a wide range 
of dark matter mass has been studied when the mass splitting takes <5 = 1,4,40,100 
GeV. We observe that the relic density changes considerably by varying the mass 
splitting <5 when moM % 60 GeV or ttidm > 400 GeV. 

In addition, the scalar split model predicts a gamma-ray excess for ttidm ~ 63 
GeV and 5 = 100 GeV which is in agreement in morphology and spectrum with 
the excess observed out of the Fermi-LAT data. To compute the gamma-ray flux 
which is produced by bremsstrahlung processes and the pion decay created from 
cascade annihilations of dark matter into SM final states, we have used the so-called 
generalized NFW halo profile for the dark matter density at high Galactic latitudes 
of the Milky Way galaxy. 
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A Annihilation Cross Sections 


In this section we present the relevant annihilation cross section formula for Feynman 
diagrams with only two particles in the final state. The annihilation cross sections 
into fermion pairs for the dark matter candidate, H 2 is 


Va,miV re \{H 2 H 2 ->■ //) = 


N c mj 
7 r 


( 1 - 


4m 2 , 3 

- — ) 2 


(Ai cos 2 9 + A 2 sin 2 9 — A 12 sin 9 cos 6)* 

(s ~ m 2 h ) 2 + mlrl 


,(30) 


and for annihilation into gauge bosons W± and Z is 


^nnVrel(H 2 H 2 ->■ W + W~,ZZ) = — 

2its 

4 ml 


(Ai cos 2 9 + \ 2 sin 2 9 — A 12 sin 9 cos 9) 


2 , 


(* - 0 2 + 


((s — + 8mjy-)(l — 


W • 


°h 
,2 \2 


4 mil 


) 2 + n(( s_ 2m 2 ) +8m w )(l — 

o Z o 


X 

• (31) 


The annihilation cross section for the process H 2 H 2 —>• /r/r involves three Feynman 
diagrams given in Fig. [2l The finall result reads 


(^ann Vre\(H 2 H 2 
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16v 4 f3 4 
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(32) 


where a = (Ai — A 2 ) sin 29 + A 12 cos 29 and j3 = Ai cos 2 9 + X 2 sin 2 9 — A 12 sin 9 cos 9. 
In the process H 2 (p\) H 2 (p 2 ) —>• h{p$) h(pn), the Mandelstam variables are s = 
{Pi +P2) 2 , t = (pi -P 3) 2 and u = (pi - pi) 2 - 
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